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ABSTRACT 

Many observations indicate that dark matter dominates the extra-galactic Universe, yet no totally dark struc- 
ture of galactic proportions has ever been convincingly identified. Previously we have suggested that VIRGOHI 
21, a 21-cm source we found in the Virgo Cluster using Jodrell Bank, was a possible dark galaxy because of 
its broad line-width (~ 200 kms -1 ) unaccompanied by any visible gravitational source to account for it. We 
have now imaged VIRGOHI 2 1 in the neutral-hydrogen line and find what could be a dark, edge-on, spinning 
disk with the mass and diameter of a typical spiral galaxy. Moreover, VIRGOHI 21 has unquestionably been 
involved in an interaction with NGC 4254, a luminous spiral with an odd one-armed morphology, but lacking 
the massive interactor normally linked with such a feature. Numerical models of NGC 4254 call for a close 
interaction ~ 10 8 years ago with a perturber of ~ 10 11 solar masses. This we take as additional evidence for 
the massive nature of VIRGOHI 21 as there does not appear to be any other viable candidate. We have also 
used the Hubble Space Telescoped to search for stars associated with the H i and find none down to an / band 
surface brightness limit of 31.1 ± 0.2 mag arcsec -2 . 

Subject headings: dark matter - galaxies: individual (VIRGOHI 21) - radio lines: galaxies 



1. INTRODUCTION 

The ability of a galaxy to form stars depends critically on 
the fraction of its mass that forms the baryonic gas disk and 
the temperature of this gas. High densities shield the gas 
from ionising radiation, allowing it to cool while higher in- 
tensity ionizing radiation keeps the gas hot. Models predict 
that galaxies can form with gas column densities that prohibit 
star formation yet provide some self shielding from the ioniz- 
ing background (Davies et al. 2006). Such dark galaxies are 
potentially detectable by blind 2 1 -cm surveys of the sky. The 
possible existence of dark galaxies has previously been dis- 
cussed by Jimenez et al. (1997), Hawkins (1997), Verde et al. 
(2002), and Davies et al. (2006); see also Taylor & Webster 
(2005) for an alternative view. 

Objects detected at 21-cm but with no optical counterparts 
have been known about for many years; these include high 
velocity clouds (Wakker & van Woerden 1997), the Leo Ring 
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(Schneider et al. 1983), and various gas clouds close to bright 
galaxies (Kilborn et al. 2000; Boyce et al 2001; Ryder et al. 
2001). However, none of these objects have the characteris- 
tics of a galaxy, i.e. detectable emission over galaxy-sized 
spatial scales and a velocity structure consistent with a rotat- 
ing and gravitationally bound disk. Recently two candidate 
dark galaxies have been reported: VIRGOHI 21 (Minchin et 
al. 2005) and HVC Complex H (Simon et al. 2006). A third 
possibility is the Hi cloud associated with the Local Group 
galaxy LGS3 (Robishaw et al. 2002). 

In this paper, we present high resolution H i observations 
of VIRGOHI 2 1 , which we believe add further support to our 
hypothesis that it is a dark galaxy. We also describe HST 
observations that were specifically designed to search for red 
giant stars at the distance of the Virgo Cluster and so place 
very faint limits on the surface density of stars that might be 
associated with VIRGOHI 21. 

2. 21-CM OBSERVATIONS AND ANALYSIS 

The new 21-cm data were taken in March 2005 at the West- 
erbork Synthesis Radio Telescope (WSRT) in two full 12- 
hour syntheses and reduced using the Miriad package. The 
data were nagged for shadowing and on two of the fourteen 
25-meter antennae one polarization was nagged due to prob- 
lems with the gain. A spectral bandwidth of 10 MHz covered 
the velocity range 930 - 3070 kms -1 . Removal of the nois- 
ier end channels left 230 useful channels of width 8.2 kms 1 
each, giving a velocity resolution of 10 km s _1 over the range 
980 to 2890 kms -1 . Continuum removal was carried out in 
the UV plane using uvlin. The standard source 3C147 was 
used for calibration. Cleaning used a robust setting of 1, close 
to normal weighting (Briggs 1995). The cleaned cube was 
Gaussian smoothed spatially and Hanning smoothed in veloc- 
ity. This was then used as a template for regions where flux 
was present in the data in order to carry out a deeper cleaning 
of the dirty image, which gave the final cube used in the anal- 
ysis. The synthesized beam was 99" x 30" in size (extended 
north-south) and the noise was 0.3 mJy beam 1 channeL 1 , 
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giving a 5cr column-density limit for sources 25 kms wide 
of 2 x 10 19 atoms cm 1 . Fig. shows a neutral hydrogen 
(H i) contour map of the field superimposed on a negative op- 
tical image. VIRGOHI 21 is the elongated structure in the 
center (which is at about 2000 kms -1 ). A faint bridge can 
be seen stretching down to the prominent spiral NGC 4254 
(2400 kms -1 ) while the other two sources, NGC 4262 (1350 
kms -1 , upper left) and the faint galaxy 'C' (1750 kms -1 , im- 
mediately to the left of VIRGOHI 21) appear unconnected. 
Fig. QJ) shows the velocity-declination projection of the full 3 
dimensional data cube. Now NGC 4254 is at the bottom right 
while VIRGOHI 2 1 is the angular structure in the center with 
'C to its left. Far more detail can be seen in the animation 
of the data cube, Fig. [2] (available in the electronic edition of 
the Astrophysical Journal, only the first frame is shown in the 
print edition). The bridge between VIRGOHI 21 and NGC 
4254 is clear, as is the lack of any connection between VIR- 
GOHI 21 and either NGC 4262 or galaxy 'C The apparent 
alignment of NGC 4262 with VIRGOHI 21 and NGC 4254 in 
Fig. Wp is merely a consequence of the particular projection 
shown. The bridge is ~ 25 arcmin long and stretches from 
2250 kms -1 at +14°28' at the low-velocity (western) edge of 
NGC 4254 to 1900 kms -1 at+14°41'. Here it meets a section 
of relatively high column-density H i which is flat in velocity 
until +14°46' where there is a strong velocity gradient in the 
opposite direction from that seen in the bridge, rising to 2100 
kms -1 at +14°49'. 

Fig. |3] is a blow-up of the source region superimposed 
on a far deeper CCD optical image. The CCD image was 
made by combining two B-band images: a 750 s Isaac New- 
ton Telescope (INT) Wide Field Camera (WFC) survey im- 
age taken during photometric conditions in March 2000 and 
our own 600 s image taken with the same instrument during 
non-photometric conditions in May 2004. Both sets of data 
were reduced using the Wide Field Survey pipeline 10 which 
includes de-biasing, bad pixel replacement, non-linearity cor- 
rection, flat-fielding and gain correction. The images were 
combined using the Starlink task MakeMos in the CCDPack 
package and photometry was taken from the WFC survey im- 
age, which has a photometric accuracy of 0.05 mag. Com- 
bined, the two images have a lcr sky noise of 26.7 mag 
arcsec -2 which has been improved to 27.5 mag arcsec -2 by 
binning into 1 arcsec pixels. This gives a surface brightness 
detection limit of about 27.5 mag arcsec -2 for an object of di- 
ameter 10 arcsec. We find no optical counterpart down to this 
surface brightness and size limit. 

There is a small, faint galaxy labeled 'A' superposed on the 
highest H i contour (10 20 cm -2 ) at declination +14°47.4'. An 
optical spectrum from the 6.5-m MMT in Arizona shows that 
it is at a redshift of z — 0.25 and is therefore unconnected 
with VIRGOHI 21 (see Fig. g). The 17th magnitude galaxy 
'C to the left at +14°45' is an H i point-source at this resolu- 
tion. By comparison, VIRGOHI 21 is an extended structure 
in both dimensions rather than a collection of discrete com- 
pact clouds. The velocity-declination plot (right) shows the 
complex kinematic structure of VIRGOHI 2 1 . The most re- 
markable feature is the tilted portion between +14°46' (1900 
kms -1 ) and +14°49' (2100 kms -1 ) which resembles the sig- 
nature of an edge-on rotating disk (e.g. Kregel et al. 2004). 
A tentative detection of gas further to the north (in Fig. 3b) 
could either be part of the interaction or clumps further out in 
the disk. Note the lack of connection between galaxy 'C and 

10 http://www.ast.cam.ac.uk/~wfcsur/technical/pipeline/ 



VIRGOHI 21. 

Figs. [5] and [6] show the channel maps from the WSRT ob- 
servations. Fig. shows the whole VIRGOHI 21 - NGC 
4254 system, including the bridge connecting the two which 
is only a 3<x detection in the individual channels but, taken 
over all channels, is definitely present (see Fig. [Tb and the 
animation of the datacube, Fig. Fig. |6]shows the cube be- 
tween RA 12 ft 17 m 46' s and 12%8 m 03 J over the velocity range 
of VIRGOHI 21. It can be seen that the section identified in 
Fig. [3J3 as the 'Rotating Disc' undergoes a smooth shift north- 
wards with increasing velocity, which resembles (taking into 
account the difference in beam sizes) that observed in edge-on 
rotating systems (e.g. Uson & Matthews 2003; Rupen 1991). 
Note that this section of VIRGOHI 21 is unresolved in in- 
dividual channels, leading to a broadening of the linewidths 
down each line of sight to ~ 100 kms -1 and a lowering of the 
column-density as the H 1 is smeared by the beam. 

3. HST OBSERVATIONS 

To place a better limit on the presence of stars associ- 
ated with VIRGOHI 21 we have obtained HST data with 
a total exposure time of 23142 s. The ACS Wide Field 
Camera was used with the /-band F814W filter. The point- 
ings were centered on the H i-flux-weighted position of 
12n7 m 53 s , +14°45'25" (J2000), with the ACS field of view 
matching the minimum H 1 size (3.5 arcmin) very closely. 
Each orbital exposure was cr-split in two for cosmic ray re- 
moval and was slightly offset to permit drizzling. 

Preliminary data reduction was done by the Space Tele- 
scope Science Institute with the standard pipeline. The flat- 
fielded images have been combined with the Iraf task Mul- 
tidrizzle to obtain a distortion-corrected and drizzled final 
image. The pixel size in the resulting images is 0.05 arcsec. 

With this depth and spatial resolution we should be able to 
resolve individual stars in VIRGOHI 21 down to / ^ 28 and 
thus to detect the tip of the Red Giant Branch which at the dis- 
tance of the Virgo Cluster would be close to 27 in /. The pres- 
ence of individual stars would not conflict with the ground- 
based INT observations: the B-band surface-brightness limit 
of 27.5 mag would (for B - I 2) equate to around 0.5 de- 
tectable RGB stars per square arcsecond, or 2500 - 3000 over 
the area of VIRGOHI 21. 

Stellar photometry was performed with the Daophot pack- 
age in Iraf (Stetson 1987). A preliminary selection of stars 
with a threshold of 3.5 sigma was performed with the au- 
tomatic star-finding algorithm Daofind and their magnitudes 
obtained with the task Phot using an aperture radius of 3 pix- 
els. A spatially variable model point-spread function (PSF) 
- with a full width at half maximum of 2 pixels - was built 
with the task Psf by selecting several isolated and bright stars. 
Finally, Allstar was used to fit the model PSF to the stars 
in the input list and to produce the final catalog. The conver- 
sion to Johnson / magnitudes and the aperture correction were 
done following Sirianni et al. (2005). Artificial star tests were 
used to determine the completeness of our selection; the 50% 
completeness limit is found to be / = 27.9 mag. 

We found that the photometry catalog produced with 
Daophot/Allstar (~600 objects) included bad measurements 
due to misidentifications such as parts of bright background 
galaxies or unresolved background galaxies whose contribu- 
tion to the luminosity function is significant at fainter magni- 
tudes (see Durrell et al. 2002). Given the long exposures re- 
quired by such observations we did not ask for additional time 
to look at another field to determine the background density 
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Fig. 2. — This figure is available as a color mpeg animation in the electronic 
edition of the Astrophysical Journal, only the first frame is shown here. The 
figure shows an animation of the WSRT data cube, oriented with increasing 
velocity upwards and rotating around this axis. VIRGOHI 2 1 can be seen in 
the center joined by a stream of neutral hydrogen to NGC 4254 (above and to 
the side). Object C is below VIRGOHI 21 and NGC 4262 is the ring structure 
near the bottom. The vertical column is the continuum source 1214+1456 
(White & Becker 1992) which was incompletely removed by uvlin. The 
faint vertical striping is due to incompletely removed continuum from M87, 
which is detected in the sidelobes of the array. 



Fig. 1. — (a) Hi contour map of the 21-cm observations, superimposed on a 1 square degree negative Digitized Sky Survey image. Contours are from 2.5 X 10" 
to 2 X 10 20 cirT 2 at intervals of 2.5 X 10" cnT 2 .(b) Shows the declination-velocity projection of the data cube. More detail can be seen in the online animation 
of the datacube (Fig.ff). 

relatively low numbers of stars in each magnitude bin. If we 
assume this TRGB to be real, despite the low signal to noise 
ratio (around 2cr on the Sobel-filter output), we can set the 
most stringent limit on the surface-brightness of VIRGOHI 
21. 

Assuming (very generously) that all of the fainter stars are 
also at the distance of the cluster (119 with / = 27.2 - 27.9) 
they have a combined magnitude of / = 22.3, giving a sur- 
face brightness of 33.8 / mag arcsec~ 2 . Using the formula of 
Durrell et al. (2002) gives a correction for only having the 
top 0.7 mags of the RGB of 1.5 + 0.1 mag arcsec~ 2 . We do 
not detect any Asymptotic Giant Branch stars; based on the 
work of Durrell et al. these could contribute another 0.2 + 0.1 
mag arcsec~ 2 . Overall, this gives an upper limit for the to- 
tal /-band surface brightness of 31.1+ 0.2 mag arcsec~ 2 . For 
normal galaxy colors of B - I — 1.5-2, this equates to an 
extremely faint B-band surface brightness of 32.4 - 33.3 mag 
arcsec~ 2 . 

These results can be compared with studies of Virgo Clus- 
ter inter-galactic stars (Durrell et al. 2002) and those used to 
determine the color-magnitude diagrams of dwarf galaxies be- 
longing to the cluster (Caldwell 2006). Durrell et al. imaged 
a field at about 40' northwest of M 87 («190 kpc) while the 
two fields observed by Caldwell were also relatively close to 
the center of Virgo, one at 1 .4° north («400 kpc) and the other 
at 1.1° west (^315 kpc). These fields are much closer to M 
87 at the center of the cluster than VIRGOHI 21 and the num- 
ber density of stars found is also much higher however, there 
is an apparent decrease in the number of stars as one moves 
away from the cluster center (see Fig. 9 in Caldwell 2006 and 
Fig. 6 in Durrell et al. 2002). A comparison between the 
estimated surface brightnesses of the three fields shows that 
Hi ranges from 27.7 / mag arcsec~ 2 (Durrell et al. 2002) for 
the fields closer to the cluster center to fij = 27.9 and 29.9 / 
mag arcsec~ 2 (Caldwell 2006) for the more distant ones. Our 
/ band surface brightness of 31.3, at a projected distance of 
1 . 1 Mpc from the cluster center (M87), is thus consistent with 
the known population intra-cluster stars in Virgo and so there 
is no evidence to suggest that our detections are associated 
with VIRGOHI 21 in any way. 



of star-like objects; instead, we have removed objects from 
our list by individually examining their radial profiles and re- 
jecting those with a FWHM > 2.5 - 3 pixels, corresponding to 
the value of the PSF built with the brightest stars in the field. 
The final list consists of 28 1 stars and the corresponding lu- 
minosity function is shown in Fig. [7] To estimate the con- 
tamination due to foreground stars, we compared the number 
of detections with a model of the expected number of Milky 
Way stars in the position of the Virgo Cluster (Castellani et 
al. 2002; Cignoni et al. 2003). This model shows that the 
contribution of foreground stars is only significant at brighter 
magnitudes (/ < 26.5) where the number of stars detected 
is small (Fig. |7); at fainter magnitudes the larger number of 
detections makes the foreground unimportant. 

To search for the Tip of the Red Giant Branch (TRGB) we 
have applied a Sobel filter to the binned luminosity functions. 
The luminosity function shows a slope change at / = 27.2+0. 1 
mag, which we identify as the possible tip of the red giant 
branch. This would place these stars at a distance of 16.6 - 
18.2 Mpc, which is consistent with them being in the Virgo 
Cluster (d=16 Mpc), though the statistics are poor due to the 
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Fig. 4. — Optical spectrum of object 'A' taken with the 6.5-m MMT. The H 
and K lines give a redshift of z = 0.2535. 

4. DISCUSSION 

In the previous sections we have described our new obser- 
vations of VIRGOHI 21. In this section we try and put them 
in context and particularly we want to assess the evidence for 
and against the dark galaxy hypothesis. 

1 . If attributed to gravitation, changes in velocity of galac- 
tic size over galactic scales, as seen here, require masses 
of galactic proportions. On dimensional grounds the 
abrupt change in velocity AV (~ 200 kms~') seen in 
VIRGOHI 21 over a conservative length-scale Ax (~ 14 
kpc =a 5 x 10 22 cm 3 arcmin at the Virgo Cluster dis- 
tance of 16 Mpc) implies a mass M > (AV) 2 Ajc/G 
jqIO-ii so j ar masses if this is due to an object rotat- 
ing in dynamical equilibrium. If the velocities are not 
due to rotation then changes in velocity like this can 
be observed if, for example, we are observing an arc 
of gas and our view is almost edge-on. These geo- 
metrical effects can lead to large velocity widths that 
may resemble rotation (e.g. Bournaud et al. 2004). 
The problem with this latter scenario (as also pointed 
out by Vollmer et al. (2005) is that the perturbing 
galaxy should be moving in the plane of the arc and 
so should be projected on the sky very close to it. We 
can simply quantify what we mean by this. Imag- 
ine two galaxies with radial velocities, V\ and Vj at 
either end of an approximately linear tidal bridge of 
physical length d pitched at an angle 9 to the line of 
sight. A telescope pointed towards it has a transverse 
beam diameter of b at the bridge. The only signif- 
icant gas motions within the bridge will be stream- 
ing velocities along its length. From end to end of 



the bridge the radial velocity difference is | V2 - V\\ 
while within the telescope beam the measured velocity 
width, AV, will be (b/d sin 6) x \V 2 - Vi\. But, bridges 
of any appreciable size arise only when the total ve- 
locity difference between the interacting galaxies i.e. 
V2 - Vi|/ cos6>, is of order the circular velocity V c in 
the galaxies involved (Toomre & Toomre 1972). It fol- 
lows that AV/V C - {b/d)cot8. Thus broad line widths 
AV ^ V c , as here, can only be seen within a beam if 
both interactors appear to lie within, or very close to 
the beam i.e. dsinO < bcosG < b. We can find no 
potential perturber within the Arecibo beam (b 3'.6) 
or close to it. Thus we do not believe VIRGOHI 21 can 
be tidal debris from an interaction of this nature. 

2. The existence of the H 1 bridge indicates that NGC 4254 
has undergone some kind of interaction. NGC 4254 is 
a luminous one-armed spiral galaxy sufficiently pecu- 
liar to have attracted several studies (lye et al. 1982; 
Phookun et al. 1993; Vollmer et al. 2005). According 
to Vollmer et al. 2005) the morphological peculiarities 
of NGC 4254 can be explained by an interaction with 
a perturbing mass of ~ 1O U M . As it is now clear that 
VIRGOHI 2 1 is involved in this interaction, it should be 
considered a candidate for the perturber, particularly as 
no other candidate can be easily identified. The Vollmer 
et al. simulation implies an interaction 3 x 10 8 years 
ago. As the projected length of the bridge is 120 kpc 
this would imply it has been drawn out at a projected 
speed of 390 kms" 1 , which is comparable to the radial 
velocity difference between NGC 4254 and VIRGOHI 
21 of 400 kms 1 . Thus if VIRGOHI 21 were massive 
enough (~ 10 11 M Q ), it could be the cause of the tidal 
bridge. 

3. If object C, the Hi galaxy just to the east (left) of VIR- 
GOHI 21 (Fig. O were involved its mass (see (a) above) 
should be ~ 10 n M o . However, its measured Hi ve- 
locity distribution and size suggest a probable mass 
< 1O 9 M , while it has a luminosity of only 10 8 L Q and a 
Mhi of 2 x 10 7 M Q . Galaxy 'C appears to be two orders 
of magnitude too undermassive and underluminous to 
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Fig. 5. — Channel maps of the VIRGOHI 21 -NGC 4254 system. Contours are at -1 mJy (— 3cr), 1 mjy (3cr), 2 mJy (6cr) and 3 mJy (9cr) (colored orange, red, 
magenta and blue in the electronic edition). The beam size is shown by the ellipse in the bottom-left and the velocity of each channel (in km s~') is given in the 
top-left. See the electronic edition of the Journal for a color version of this figure. 



explain the one-armed spiral structure of NGC 4254 
and the H i bridge connecting to it. So, it seems very 
unlikely that 'C has sufficient mass to cause a pertur- 
bation of over 200 km s -1 to the tidal stream, and VIR- 
GOHI 21 passes to its west, not between it and NGC 
4262 as might be expected if it had pulled the stream 
westward. The radial velocity of 'C means that it must 
be moving past the stream at a velocity (relative to the 
putative undisturbed velocity of the stream) of at least 
350 kms~'. At this speed, it would not have stayed 
close to the stream long enough to have severely dis- 
turbed it. It would also be expected that if 'C were in- 
volved in the interaction then there would be gas falling 
onto it, but this does not appear to be the case as there 
is no gas seen between VIRGOHI 21 and C. The Hi 
in 'C is unresolved, implying that it is confined to the 
area of the optical galaxy; it shows no signs of infalling 



gas streams or any other disturbance. Galaxy 'C is in 
our a minds a possible, but highly unlikely participant 
in this interaction. 

4. Recently Bekki et al. (2005) have published numerical 
simulations on the basis of which they argue that "VIR- 
GOHI 21 ... is likely to be tidal debris rather than a 
'dark galaxy' ". They model high-speed, hyperbolic in- 
teractions within cluster potentials (at 350 and 700 kpc 
from the cluster center) over 2 Gyr in order to draw out 
very long tidal tails where the H i and stars have become 
separated so that the average surface brightness within a 
100 kpc 2 box is > 30B mag arcsec^ 2 (L B < 6.9xl0 8 L o ), 
assuming M+/Lb = 4. The Bekki et al. stellar den- 
sity can be converted to /-band if we assume a value 
for Mi, I Li. Estimates of this vary (e.g. Portinari et al. 
2004; Vallejo et al. 2002; McGaugh & de Blok 1997) 
but are generally close to 1 and always less than 2. 
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Fig. 6. — Channel maps of the central region of VIRGOHI 21, contours as in Fig. H] The cross indicates the peak of the zeroth moment map, which is at the 
center of the rotating component (see Fig. [3). It can clearly be seen that this component (declination range 14° 46' to 14° 49') starts south of this point and moves 
northwards across it with increasing velocity (as might be expected for rotation). See the electronic edition of the Journal for a color version of this figure. 




Fig. 7. — The luminosity function of the final list of stars. The dotted line 
shows the output of the edge detection filter used to detect the tip of the RGB. 
The maximum of the function indicates the position of the tip. The vertical 
line at / = 27.9 mag indicates the 50% completeness limit. 



We therefore adopt a very conservative (in the sense of 



producing the lowest /-band surface-brightness) value 
of M-t/Li = 2. This gives a value of B - I for the 
surface-brightness estimates given in Bekki et al. of 
0.75. For the two double-tailed models presented (Ml 
and M2) this gives average /-band surface-brightnesses 
of 32.7 and 33.0 mag arcsec~ 2 , while for the two single- 
tailed models presented (M3 and M4) this gives aver- 
age /-band surface-brightnesses of 29.5 and 29.9 mag 
arcsec~ 2 . The two-tailed models are inconsistent with 
our Hi data and the Hi data of Phookun et al., which 
clearly show a single tail from NGC 4254, while the 
single-tailed models give considerably higher surface- 
brightnesses (1.2 - 1.6 mags, or 6 - 8 <x) than the upper 
limit from our HST observations. In addition to this, 
Bekki et al. find that the velocity fields produced in the 
dark clouds by their models do not resemble rotation, 
concluding that "Velocity fields are the key observa- 
tional tools to help us determine whether gas clouds are 
(unbound) tidal debris or are self-gravitating systems 



21 -cm observations of VIRGOHI 21 



7 



embedded within a massive dark matter halo." From 
both the surface-density of stars expected in a single- 
tail interaction and from the rotation curve, it appears 
that high-speed, hyperbolic interactions within a cluster 
potential like those proposed by Bekki et al. are unable 
to reproduce our new observations of VIRGOHI 21. 

5. NGC 4262, the spiral to the North-East, does not ap- 
pear to be involved. This is not surprising because it is 
in totally the wrong place to give rise to a tidal feature 
like that observed (NGC 4262, NGC 4254 and VIR- 
GOHI 2 1 form a well separted triangle on the sky, see 
also Vollmer et al. 2005). Also, the radial velocity be- 
tween NGC 4254 and NGC 4262 is far too large (900 
kms~') to generate tidal features such as bridges and 
tails. Furthermore, with such a large radial velocity dif- 
ference between NGC 4254 and NGC 4262, any inter- 
action must be at high speed and, as seen above from 
the simulations of Bekki et al., such interactions are 
not capable of reproducing our observations as even the 
most isolated gas clouds formed in single-tail systems 
have surface-brightnesses 1.5-2 magnitudes brighter 
than the limit set by our HST observations. NGC 4262 
is a SO galaxy with an H i ring that is inclined with re- 
spect to the optical disk (Krumm et al. 1985). Rrumm 
et al. present higher-resolution Hi observations than 
our data and they say that if the H i has come from an 
interaction, "the narrowness of the ring and the undis- 
turbed appearance of the optical disk suggest that the 
donor was a gas-rich dwarf galaxy or intergalactic H i 
cloud, rather than a massive spiral". Thus we find it 
very difficult to see how NGC 4262 can have given rise 
to the features we have now observed. 

6. The components of VIRGOHI 21 are connected both 
spatially and in velocity, making it exceedingly improb- 
able that they could be chance superpositions of smaller 
hydrogen clouds, while the bridge to NGC 4254 is not 
explained by this hypothesis. 

7. Phookun et al. (1993) suggested that the distortion of 
NGC 4254 could be due to infalling gas-clouds. Could 
VIRGOHI 21 be a tail, similar to those seen in UGC 
10214 (Briggs et al. 2001) or NGC 4038/9 (Gordon 
et al. 2001)? In NGC 4038/9, in particular, there is 
what appears to be a tidal dwarf near a strong concen- 
tration of H i at the end of the southern tail, around 60 
kpc from the center of the system (for a distance of 13.8 
Mpc; Saviane et al. 2004). However, the major prob- 
lem with this hypothesis is that NGC 4254 apparently 
shows no sign of having undergone a recent merger. It 
would be expected that a merger violent enough to have 
thrown out an arm over 100 kpc in length would leave 
some imprint on the optical disk, but this is not the case 
here. Phookun et al. found that "The galaxy is under- 
going an interaction that is not so strong or violent as 
to disrupt the disk." While this fits with the hypothesis 
of VIRGOHI 21 as a dark galaxy, playing a role similar 
to that of the companions of NGC 4027 or NGC 4654 
in exciting the m - 1 mode, it does not fit the merger 
scenario. 

8. Oosterloo & van Gorkom (2005) argue that another H i 
cloud in the Virgo Cluster, VIRGOHI 4 (Davies et al. 
2004), is caused by ram-pressure stripping from NGC 



TABLE 1 
Properties of the putative dark disk 



14 kpc 
100 km s -1 
4 X 10 8 years 
2 X 10 10 M o 
2 X 10~ 2 g cm" 



Diameter 2R, [from +14° 46' to 49'] 
Circular Velocity V c , [(2100 - 1900)/2] 
Spin Period P, [2nR/V c ] 
Total Mass M T , [RV 2 /G] 
Face-on Mass-density [Mr/nR 2 ] 
Hydrogen Mass Mm, [F H i = 0.6 Jykms" 1 ] 3 X 10 7 M o 
Face-on gas density Nfji, [Mm /nnifjR 2 ] 3 X 10" cirT 2 

Luminosity limit Lg, [jig > 32.8 mag arcsec~ 2 ] < 1.9 X 10 5 L o 
Total Mass to Blue Light Ratio M dy ,JL B > 10 6 M /L 

[assuming a disk 0.5 arcmin X 3 arcmin] 



TABLE 2 

H I FLUXES AND MASSES OF THE DETECTED COMPONENTS 



F HI (Jy kms ') 



NGC 4254 


91.5 


5.5 x 10 9 


NGC 4262 


7.7 


4.6 x 10 8 


Object C 


0.3 


2x 10 7 


'Rotating Disc' 


0.6 


3 x 10 7 


(14° 46' < d < 14° 49') 






Low col. density bridge 
(6 < 14°41') 


1.3 


8 x 10 7 






High col. density bridge 


0.8 


5 x 10 7 


(14°41' < 6 < 14° 46') 







4388 due to an interaction with the hot-gas halo of 
the M86 sub-group and suggest a similar origin could 
be possible for VIRGOHI 21. This would explain the 
bridge without the need to invoke a second galaxy, ei- 
ther interacting or merging. However, we cannot see 
how this can give the steep, reversed velocity gradient 
seen in VIRGOHI 21, nor does it explain the distortion 
to the optical disk of NGC 4254. 

9. Single-armed spirals are normally the result of inter- 
actions with close-by massive companions (lye et al. 
1982; Phookun et al. 1993); the lack of any visi- 
ble companion to NGC 4254 has thus led to a num- 
ber of observations and dynamical models. As men- 
tioned above, recent numerical models by Vollmer et 
al. (2005) suggest that NGC 4254 "had a close and 
rapid encounter with a 10 n M o galaxy ~ 250 Myr ago. 
The tidal interaction caused the spiral structure. . . ". 
Phookun et al. (1993), in a VLA study of the galaxy, 
find a trail of gas leading away from it (their Figure 5) 
which can be identified from its direction and velocity- 
gradient as the NGC 4254 end of the 'bridge' between 
the two objects. This supports our inference that VIR- 
GOHI 21 is the aforesaid ~ 10 n M o mass which caused 
the peculiarities in NGC 4254 - particularly so given 
the lack of any other candidate. 

If our hypothesis is correct that this is a dark, gravitationally 
bound, edge-on rotating disk then its properties are as pre- 
sented in Table Q] Judging from visible disk galaxies, whose 
masses continue to rise beyond their H i edges (e.g. Salucci & 
Persic 1997), the full size and mass of such a disk could eas- 
ily reach the ~ lO n M required by the Vollmer et al. (2005) 
simulations. The very low surface-brightness limits give an 
upper limit to the luminosity of the disk of 1.9 x 10 5 L Q , giv- 
ing a lower limit for M</ v „/L# of at least 1O 6 M /L , whereas 
normal galaxies have values < 50. 

It is clear from our new observations that VIRGOHI 21 is 
intimately linked to NGC 4254. The important question is 
whether it is the cause of the HI bridge or the effect of some 
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previous encounter. We believe the weight of the evidence 
presented above lies on the side of a dark galaxy. Dark halos 
have, after all, been predicted by galaxy evolution simulations 
such as Jimenez et al. (1997), Verde et al. (2002) and Davies 
et al. (2006). That none have been confirmed so far is most 
probably due to instrumental limitations. VIRGOHI 21 itself 
was missed in the normal HIPASS-sensitivity observations of 
this area, and only picked up with ten times longer integra- 
tions in the VIRGOHI survey. In any case single-dish angular 
resolution makes it all too easy to confuse a dark cloud with a 
bright object clustered with it (Davies et al. 2006). The recent 
optical identification program for HIPASS sources (Doyle et 
al. 2005) can rule out only 'isolated', i.e. unclustered, dark 
galaxies. 

If we accept the total mass given in Table 1 then an odd- 
ity about VIRGOHI 21 is its shortage of baryons compared 
to visible galaxies. Visible disk galaxies typically have, in 
the form of stars and gas, ~ 10 percent of their dynamical 
mass in baryons, a figure consistent with cosmological mod- 
els which predict ~ 7 times as much Dark Matter as bary- 
onic in the Universe at large (Salucci & Persic 1997). Had 
VIRGOHI 21 ten percent of its mass in baryons it would con- 
tain ~ 2 x 10 9 solar masses which would, if it were once in 
the form of H i spread evenly across the existing dynamical 
disk, yield a column density ~ 10 21 atoms cirT 2 - sufficiently 
above the "star-formation threshold" of ~ 10 20 cirT 2 (Martin 
& Kennicutt 2001), or above the Toomre criterion for gravi- 
tational instability in a rotating disk (Toomre 1964), to form 
stars (Jimenez et al. 1997; Verde et al. 2002). Since it appar- 
ently contains no stars it must have a minimum of ten times 
less baryons - a phenomenon often apparent in dwarf galaxies 
(Mateo 1998). Thus 90 percent of the baryons one might have 
expected to find in a typical galaxy appear to be missing and 
so it is not difficult to see how a galaxy like this may remain 
dark. 



5. SUMMARY 

In this paper we have presented new Westerbork high reso- 
lution 21 -cm observations of a HI source (VIRGOHI 21) that 
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we have previously observed at Jodrell Bank and Arecibo. 
These observations clearly show that VIRGOHI 21 has played 
a part in some form of interaction with the bright spiral galaxy 
NGC 4254. Whether VIRGOHI 21 is the result of some pre- 
vious interaction with a third object or it is the interacting 
object itself is still open to conjecture. We have argued in this 
paper that the weight of evidence now resides on the side of 
the latter rather than the former conjecture. This is primar- 
ily because the new observations show a long drawn out tidal 
feature which has an abrupt velocity change and we cannot 
find a candidate galaxy to have caused it. Tidal interactions 
between optically bright galaxies often draw out stars as well 
as gas. We have also used HST to search for stars associated 
with VIRGOHI 21 and find none down to very faint surface 
brightness levels. We look forward to new observations and 
numerical models that will add further weight to either side of 
this debate. 
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